Entanglement Entropy in an Antiferromagnetic Heisenberg Spin 

Chain with Boundary Impurities 

Jie Ren 1 *, Shiqun Zhu 1 " 1 ", and Xiang Hao 2 
1 School of Physical Science and Technology, Suzhou University, 
Suzhou, Jiangsu 215006, People's Republic of China and 
2 Department of Physics, Suzhou University of Science and Technology, 
Suzhou, Jiangsu 215011, People's Republic of China 

Abstract 

The effects of boundary impurities on the entanglement entropy in an antiferromagnetic Heisen- 
berg opened spin-1/2 chain are investigated. The method of density-matrix renormalization-group 
is used to obtain the bipartite entanglement. The entropy increases when the length of the sub- 
system increases. It will approach to a constant when system length is very large. With the same 
impurity interaction, qutrit impurities of spin-1 can increase the entanglement entropy. 
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I. INTRODUCTION 



Recently, the entanglement has been recognized as an important resource of some quan- 
tum mechanical phenomena, such as quantum teleportation, quantum cryptography, quan- 
tum computation, and violation of Bell's inequality p~H3]. Many investigations show that 
entanglement exists naturally in the spin chain when the temperature of system is at zero. A 
useful many-body entanglement measure of a pure state is the von Neumann entanglement 
entropy [3], which can also quantify quantum phase transition [oTfTT]. The bipartite entangle- 
ment in systems of atomic Bose-Einstein condensate was studied [T2l IT3] . The entanglement 
entropy in the antiferromagnetic Heisenberg XX chain and Ising model was investigated 
[13]. In isotropic antiferromagnetic Heisenberg model, the universal form of the entropy is 
predicted by [T5] 

c 

Sl = 7^log 2 L + k, (1) 

where c is the central charge and A; is a non-universal constant. For a spin chain of open 
boundary condition, the analogous formula c/3 should be replaced by c/6 for a part of length 
L in an infinite one-dimension system [161 HZ] ■ The analogous formula is dependent on the 
boundary conditions of the block and the rest of the chain. If the block has two boundaries 
with the rest of the chain then the factor is c/3, while if the block has just one boundary as 
in the case of a block consisting of the first adjacent spins of a semi-infinite chain, then the 
factor is c/6. Recently, it is shown that a feeble central bound defect JTS] or single impurity 
in the boundary [19] has strong influence on the entropy, though the entropy measures the 
mutual coupling of the two parts of a system in wave function. A weak transverse boundary 
magnetic field impurity [20] and domain walls [21] generated by antiparallel magnetic field 
have different effects on the entanglement entropy. Moreover, impurities show their strong 
influences on spin correlation function [2"2"l 123] . It would be interesting to investigate the 
effects of boundary impurities on the entanglement entropy in an opened antiferromagnetic 
Heisenberg spin chain. 

In this paper, the entanglement entropy of a spin— 1/2 antiferromagnetic Heisenberg chain 
with boundary impurities located at two ends is investigated. In scetion II, the Hamiltonian 
of an antiferromagnetic Heisenberg spin— 1/2 chain is presented. By using the method of 
density-matrix renormalization-group (DMRG) [231 125] , the entropy of ground state is cal- 
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culated and the effects of the impurities are analyzed in section III. A discussion concludes 
the paper. 



II. HAMILTONIAN OF HEISENBERG SPIN CHAIN 

The Hamiltonian of a spin— 1/2 Heisenberg opened chain with boundary impurities at two 
ends can be written as 

N-2 

H=J2 J ^+i + a(Sj 2 + S N -J N ). (2) 

i=2 

Where the coupling exchange J > corresponds to the antiferromagnetic case, Sj are spin 
operators, N is the length of the spin chain. The coupling exchange a is impurity interaction. 
For simplicity, J = 1 is assumed in this paper. 

The entropy is used as a measure of the bipartite entanglement. If \Gs) is the ground 
state of a chain of N qubits, a reduced density matrix of L contiguous qubits can be written 
as 

p L = Tr N - L \Gs)(Gs\. (3) 

The bipartite entanglement between the right-hand L contiguous qubits and the rest of the 
system can be measured by the entropy 

S L = -Tr{p L \og 2 p L ). (4) 

One of the properties of the entropy of a block of the system can be given by 

Sl = Sn-l, (5) 
since the spectrum of the reduced density matrix p^ is the same as that of Pn-l- 

III. ENTANGLEMENT OF ENTROPY WITH IMPURITIES 

In order to calculate the entropy accurately using the method of DMRG, the length of the 
spin chain needs to be relatively long. The length of the spin chain is chosen to be iV = 256. 
The total number of the density matrix eigenstates held in the system block is m = 128 in 
the basis truncation procedure. 
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To check the accuracy of the results from the method of DMRG, the open boundary 
condition without impurities of a = 1 is considered. The corresponding results of finite spin 
chain, which is predicted by conformal field theory (CFT), can be considered as a benchmark. 
It can be written as 

S L = %g 2 [-sm(^L)] + A J (6) 

6 7T iV 

where c is the central charge and A is a non-universal constant [2"6"j |2"T] . There are large 
oscillations between even and odd L- value entropy. To avoid these relatively large oscillations, 
the even- value entropy is chosen. The entropy Sl between contiguous L qubits and the remain 
N — L qubits is plotted as a function of the subsystem L in Fig. 1 when N = 160, 200 and 
256. For L < 8, the results of DMRG are slightly lower than that of CFT. For L > 8, almost 
perfect agreement between the two results is obtained. For different values of N with large L, 
Sl is small for small values of N. For small L, there is almost no difference between different 
values of N. It seems that Sl approaching a constant for very large L is mainly due to finite 
size effect. The entropy Sl is also plotted as a function of log 2 [^ sin(-^L)] in the inset of Fig. 
1. It is shown that the entropy appears as a straight line whose slope is very close to c/6. 

The entanglement entropy Sl is plotted as a function of the subsystem length L for 
different values of the impurity interaction a in Fig. 2(a). It is seen that the entropy Sl 
increases with the subsystem length L and then approaches a constant when L is very large 
for a = 0.1,2.0. When a = 0.3,0.5, the entropy Sl decreases slightly, then increases and 
approaches a constant for very large L. The minimal value is 1.53 at L = 6 for a = 0.3 
and 1.16 at L = 4 for a = 0.5. When a = 0.1, Sl approaches the value about 2.4 for very 
large L. While for a = 0.3, 0.5, 2.0, Sl approaches the value about 1.6 for very large L. The 
influence of the impurity at two ends of the Heisenberg spin— 1/2 chain depends on the value 
of the impurity interaction a. For a = «o = 0.235, the strength alternation of even bond 
and odd bond in the center of the spin chain is minimized to almost close to zero [24J. For 
a < QiO) the even sublattice is favored. This induces larger value of Sl- While for a > cto, 
the odd sublattice is favored. This induces smaller value of Sl- If the value of a is less 
than ao = 0.235, the effects of the impurity on the entanglement entropy Sl is stronger. 
Therefore, the entanglement entropy Sl of a — 0.1 is much larger than the Sl of a = 0.3, 0.5 
and 2.0 [IH1 ESI IMJ- The entropy Sl is also plotted as a function of log 2 [^ sin(-^L)] in the 
inset of Fig. 2(a). It is seen that Sl is almost a straight as a function of log 2 [^- sin(-^L)] for 
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very large L. 

If Sia is the entropy with impurities at two ends and Slo is the entropy without impurities, 
the difference of the entropy ASl can be defined as 



The entropy difference ASl may also be called "impurity entanglement entropy" that is 
induced by adding impurities at two ends of the spin chain [19]. The entropy difference 
ASl is plotted as a function of the subsystem length L for different values of the impurity 
interaction a in Fig. 2(b). It is shown that ASl decreases and then approaches a constant 
when the subsystem L increases for a = 0.1,0.3,0.5. The value of ASl increases and then 
approaches a constant with increase of the subsystem length L when a = 2.0. It seems that 
the effect of impurity decreases with the increase of the subsystem L when a < ao = 1.0. 
However, the effect of the impurity increases when a > ao = 1.0. The entropy difference 
ASl is also plotted as a function of log 2 [— sin(^L)] in the inset of Fig. 2(b). It is seen 
that ASl is almost a straight as a function of log 2 [^ sin(-^L)] for very large L. Similar to 
that shown in Fig. 2(a), the entropy difference ASl of a = 0.1 is much larger than ASl of 
a = 0.3, 0.5 and 2.0. This is mainly due to the fact that small value of a < ao = 0.235 can 
induce stronger effect of impurity on ASl since the even sublattice is favored [T8l [T9l ¥IQ . 

From Fig. 2, it is seen that both values of Sl and ASl are decreases when a increases 
especially for small number of L. There are large differences of Sl and ASl for different 
impurity interactions a when L is small. If L is quite large, Sl and ASl approach to 
constants. If a = 0.3, 0.5 and 2.0, Sl approaches to about 1.6 while ASl approaches to 
about zero for quite large L. It seems that the effect of the impurity at two ends is very 
small for large L. If a = 0.1, both Sl and ASl are quite large. It seems that the small value 
of impurity interaction a < ao = 0.235 can induce strong effect on and ASl- 

The central charge c in Eq. (6) plays an important role in the measurement of the 
entanglement entropy. The central charge c can be calculated numerically by [T?t [T9"] 



ASl — Slo — Slo- 



(7) 



c(L) = 6[ 



T(L + 2) -T(L-2) J ' 



L+2 — iJL-2 



where 




(9) 
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The central charge labeled by c(L) is plotted in Fig. 3( function of the subsystem 

length L for different values of the impurity interaction a. When a = 0.1, the central charge 
c(L) increases to a peak and then decreases slowly with the increase of the subsystem length 
L. The central charge c(L) decreases and then approaches a constant with the increase of 
the subsystem L when a = 2.0. When a = 0.3, 0.5, the central charge c(L) increases and 
almost approaches a constant with the increase of the subsystem length L. The central 
charges of c(L = 6) and c(L = 4) are negative when a = 0.3 and 0.5 respectively. This 
corresponds to the minimum values of Sl shown in Fig. 2(a). For a > 0.235, the central 
charge c(L) increases with increasing value of a. For a = 0.1 < 0.235, c(L) of a = 0.1 is 
larger than that of a = 0.3, it decreases and finally approaches to that of a = 0.3 when 
L is very large. The value of c(L) of a = 0.1 is larger than that of a = 0.5 if L < 20. 
If L > 20, c(L) of a = 0.1 is smaller than that of a = 0.5. This is mainly due to the 
stronger singlet bonds on the even numbered links of the chain for a < 0.235 [23]. Since 
the central charge may clarify the behavior of the entropy for large values of subsystem, the 
central charge c(L = 80) is plotted as a function of impurity interaction a in Fig. 2(b) when 
1 L{= 80) < N/2(= 128). It is seen that the central charge c reaches a minimum value 
when impurity interaction a = 0.235. When impurity interaction a < 0.235, the central 
charge c decreases with impurity interaction a increases. The central charge c increases 
with impurity interaction a increases when a > 0.235. It approaches to 1.0 when impurity 
interaction a is close to 2.0. 

— * 

If the impurities are qutrit with spin-1 operators S', the effects of qutrit-impurities on 
the entropy can also be investigated. The entanglement entropy Sl and the difference of the 
entropy ASl are plotted in Figs. 4(a) and 4(b) respectively as a function of the subsystem 
length L for impurity interaction a and different impurities. The entropy Sl and the dif- 
ference of the entropy ASl are also plotted as a function of log 2 [^ sin(-^L)] in the insets of 
Figs. 4(a) and 4(b). Similar to that shown in Fig. 2, both Sl and ASl of spin-1 decreases 
with the increase of a. The entropy Sl increases and then almost approaches a constant 
when the the subsystem length L increases. The entropy of impurity with spin-1 is much 
larger than that with spin-1/2. For the impurity of spin-1/2, the entropies of a = 0.5 and 
a = 2.0 are almost not distinguishable when the subsystem length L is very large. While for 
the impurity of spin-1, the difference of the entropies of a = 0.5 and a = 2.0 are quite large 
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and approach a constant when L is very large. It is clear that the effects of qutrit impurity 
on the entanglement entropy are much stronger than that of qubit impurity. It seems that 
it is more easily to control the entropy of the system using qutrit impurity. 

IV. DISCUSSION 

It is clear that the impurity interaction and the impurity spin have a strong influence on 
the entanglement of the two subsystems [191 EH EH] . For pairwise entanglement between the 
impurity spin and the spin chain, the two boundary spins will have a strong tendency to 
form a singlet pair when the impurity interaction is large. This will reduce the entanglement 
between the boundary of the two spin subsystems and the rest of the system. The value 
of entanglement entropy is mainly determined by the density-matrix spectra, extremely by 
the few largest eigenvalues of the reduced density matrix [TBI EQl EE]. For qubit impurities, 
impurities can affect the entropy between two subsystems by changing the distribution of 
the reduced density-matrix spectra. If the impurities are qutrits with the same impurity 
interaction, not only the distribution of the reduced density-matrix spectra is changed, but 
also the degree of the freedom of density-matrix spectra of the subsystem is enlarged in the 
Hilbert space. This is similar to the result of the entropy with the increase of subsystems 

H5U2DH2B- 

The effects of boundary impurities on the bipartite entanglement in an antiferromag- 
netic Heisenberg open spin chain are discussed. Using the method of density-matrix 
renormalization-group, entanglement entropy is calculated for the even number subsystem. 
The entanglement entropy decreases with the increase of the impurity interaction while it 
increases with the increase of the subsystem length. When the system length is very large, 
the entanglement entropy approaches a constant due to finite size effect. The influences of 
boundary impurities with qutrit of spin-1 are much stronger than that of qubit of spin-1/2. 
With the same impurity interaction, qutrit impurities can increase the entanglement. All the 
results are dependent with the selection of even subsystem. This is shown that the entropy 
of the system with qutrit impurity can be more easily controlled. 
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FIG. 1: The entropy Sl between contiguous L qubits and the remain N — L qubits is plotted as a 
function of the subsystem length L for N = 160, 200, 256 (from bottom to top). The dashed line is 
obtained from CFT and the symbol Q is obtained from DMRG. 
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FIG. 2: The entanglement entropy Sl and the entropy difference ASl are plotted as a function 
of the subsystem length L for different values of the impurity interaction a. (a). The entropy Sl- 
(b). The entropy difference AS L . The symbols are for a = 0.1(B), 0.3(o),0.5(»), 2.0(D). 
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FIG. 3: (a). The central charge c(L) is plotted as a function of the subsystem length L for different 
values of the impurity interaction a. The symbols are for a = 0.1(B), 0.3(o), 0.5(»), 2.0(D). (b). 
The central charge c(80) is plotted as a function of impurity interaction a. 
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FIG. 4: The entanglement entropy Sl and the entropy difference ASl are plotted as a function of 
the subsystem length L for different values of the impurity interaction a and the spins 1/2 and 1. 
(a). The entropy Sl. (b). The entropy difference ASl. The symbols are for qubit impurity with 
a = 0.5(»), 2.0(o), qutrit impurity with a = 0.5(B), 2.0(D). 
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